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The assumptions underlying models based on discrete 
time usually consider all susceptible and infested in- 
dividuals to be mixed together homogeneously. This 
situation is most nearly represented by small groups of 
trees but does not hold for large stands. However, in- 
cubation and latent periods are not variable, and the in- 
festing of a tree can be considered as a relatively short 
period of time. As this model was refined, habitat type 
and volume yield factors were included. These factors 
govern tree and stand susceptibility and affect the life 
processes of beetle populations. 

One important problem with the chain binomial model 
is that substantial departures from the assumptions of 
constant incubation and latent periods and a very short 
infectious period would invalidate the model. Another 
problem is failure to properly identify the links of the 
chain. However, if a highly variable incubation period oc- 
curs, or the symptoms cannot be identified correctly, 
there is still an alternative—to base our analysis on the 
total number of cases occurring during the course of the 
epidemic. Some precision is lost when the parameters are 
estimated. However, if the number infested is large, fre- 
quencies based on this number can be calculated directly 
and will probably be more accurate than those derived 
from the probabilities of the individual chains. 


THE RATE OF LOSS MODEL 


If p is the probability of a tree becoming infested in a 
given time interval, then q = 1 — p is the probability of 
a tree not becoming infested. The probability of a tree 
becoming infested depends on the susceptibility or 
resistance of the tree, the infestivity of the beetle, the 
length of attack period, and the size of the attacking 
beetle population, as well as the environmental condi- 
tions of the stand. 


If D, is the number of trees infested at time t, then 
q?t is the probability that a specified tree will not be in- 
fested, and 1 -qP: is the probability that the tree will be 
infested. If there are G, green trees capable of being in- 
fested in the population at time t, the expected number 


of infested trees produced at the time t + 1 is G, times 

the probability of at least one tree being infested. Or: 
D,+1 = G, (1 - q?t)and G,,, = G,q>t. 

This equation provides a method of stepwise calculation 

of trees infested at successive time periods as shown in 

table 1. 

If G, = 0, all the trees are dead—no more susceptible 
trees are left—and the epidemic subsides due to food 
depletion. If D, = 0, there are no more trees successfully 
producing beetles—and the epidemic subsides. 

The Greenwood model postulates that the probability 
of a susceptible tree being infested is a constant and is 
not related to the number of infested trees. In other 
words, a susceptible tree in a stand with one infested 
tree is as likely to be attacked as the same tree sur- 
rounded by many infested trees. This is obviously not 
the case. Thus, we adopted the Reed-Frost model for 
susceptibility because it accounts for the increase in in- 
festation pressure due to the number of infested trees. 
In the Reed-Frost model, the probability of a tree not 
being infested from only one source is taken to be a con- 
stant, q. The probability of not being infested from two 
sources is thus (q) (q), and consequently from n sources 
it is qr. 

The value of q, the probability of a tree not being in- 
fested from one source, can be calculated by solving the 
equation of G,,, for q. This yields: 

q = (G,,,/G,)(1/D)) 

Theoretically, q will be a constant, but the real world 
is never constant. Thus the q for time t (q,) varies slight- 
ly with t, and may be determined for each time interval. 
However, we found a closer prediction of D,,, was ob- 
tained if several values for q, were calculated, and q was 
estimated by q, for several stands. We also found that 
precision of prediction increased with decreasing size of 
diameter classes. Estimates of tree mortality over time 
approximated true losses more closely when predicted by 
2-inch (5.1-cm) diameter classes than by larger diameter 
classes. 


Table 1.—Calculation of a theoretical epidemic from the Reed-Frost 


model (p = 0.5) 


Number Number of 
Time of dead susceptible 
period trees trees 
0 1 100 
1 5 95 
2 21 74 
3 49 25 
4 23 2 
5 1 1 
6 0 1 


Calculation of D,,, and G;., 


D, = 100 (1 — 0.95) = 5.00 = 5 
= 100 —5= 95 
D> = 95 (1 — 0.955) = 21.49 = 21 
G, = 95 — 21 = 74 
74 (1 — 0.9521) = 48.80 = 29 
G, = "74 = 49 = 25 
D, = 25 (1 — 0.9549) = 22.97 
= 25 — 23 =2 , 

= 2(1 — 0.9523) = 139 = 1 
Gs =2-1=1 
Ds = 1(1 — 0.951) = 0.05 =0 
=1-0=1 
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For high q values, peak mortality tended to be 
overestimated. By the third year, q usually becomes 
small due to the ‘“‘high-grading”’ action of the beetle in 
thinning a stand, resulting in greater overestimation of 
tree mortality. However, the critical time during an in- 
festation by the mountain pine beetle is at the point of 
change from endemic to epidemic. The value q applied to 
the larger diameter trees forecasts the pending infesta- 
tion adequately in spite of the tendency toward 
overestimation. 

The model assumes optimum conditions for the life of 
the epidemic. However, actual field conditions can cause 
beetle populations to deviate from predictions. 
Overestimation of tree mortality is not considered 
serious in most cases, particularly in the larger diameter 
classes. Epidemics usually begin in larger diameter trees 
preferred by the mountain pine beetle, and the rate of 
tree loss within these classes is critical. Thus, any factor 
that affects brood survival (such as thick phloem [food 
supply], which is correlated with larger diameters) will 
affect the rate of tree loss and, in turn, successive 
generations. 

Dispersion of the beetle is also affected by stand 
characteristics such as species, age, stocking levels, 
growth rates, and diameter class distribution; and by 
site characteristics, including habitat type, soils, eleva- 
tion, slope, and aspect. During the past decade the 
system of environmental classification by habitat type 


developed by Daubenmire for the Northern Rocky Moun- 
tain Forest Ecosystem has gained increasing acceptance 
in other areas of the West. This concept stresses use of 
the entire climax plant community as an environmental 
indicator that permits identification of environments 
(habitats) with similar biotic potentials. All en- 
vironments (habitats) with the potential to support ap- 
proximately the same mix of stable (climax) plant 
species are considered to be within the same habitat 
type regardless of successional status of the vegetation. 

Recent data from the Forest Service Northern Region 
show that the extent of lodgepole pine mortality caused 
by the mountain pine beetle varies by habitat type 
group, and by habitat type within groups. This type of 
classification provides the framework essential for 
organizing information to select alternative management 
activities for habitat types. 


MODEL TESTS AND REFINEMENT 


Data from a mountain pine beetle infestation in the 
Bechler River Drainage of Yellowstone National Park 
(Klein and others 1978) were used to predict tree loss by 
2-inch (5.1-cm) diameter classes (situation A, table 2; fig. 
1). Trees were also grouped by 6- to 12-inch (15.2- to 
30.5-cm), 14- to 16-inch (35.6- to 40.6-cm), greater than 
16-inch (40.6-cm) diameters, and total stand (table 3; fig. 
2 and 3). 


Table 2.—Predicted versus observed tree loss by year based on g;, the average probability of tree loss by 2-inch (5.1-cm) tree 
diameter class (Situation A, Klein and others 1978) 


el 


Number of trees 


Diameter Year of per acre Predicted 
class infestation nee eer I) OO 1/D, at tree loss 
6-inch 0 79.8 0.3 3.333 0.9875 0.14 

(15.2-cm) 1 79.5 0 0 1.0000 0 

2 79.5 3 3:333 .9857 14 
3 79.2 2.1 .476 .9873 99 
4 77.1 0 0 1.0000 0 
5 77.1 0 0 1.0000 0 
6 77.1 0 0 0 
Total loss: 2.7 Average: .9937 U7 
(6.75/ha) (3.175/ha) 
8-inch 0 62.7 8 1.250 0.984 0.90 
(20.3-cm) 1 61.9 8 1.250 .984 89 
2 61.1 2.7 .370 983 2.92 
3 58.4 8.1 120 .982 7.99 
4 50.3 ath 1.430 .980 63 
5 49.6 5) 2.000 -980 .63 
6 49.1 0 0 0 
Total loss: 13.6 Average: .982 13.96 
(34.0/ha) (34.9/ha) 
10-inch 0 38.8 0.8 1.250 , 0.974 1.09 
(25.4-cm) 1 38.0 11 .909 974 1.46 
2 36.9 3.9 .256 972 4.79 
3 33.0 10.6 .0943 .964 10.38 
4 22.4 1.4 714 955 1.09 
5 21.0 6 1.667 .953 44 
6 20.4 2 5.000 14 
7 20.2 0 0 0 
Total loss: 18.6 Average: -965 19.39 
(46.5/ha) (48.475/ha) 
12-inch 0 17.0 0:6 1.667 0.942 0.95 
(30.5-cm) 1 16.4 iS} .769 .938 1.91 
2 (5s 2.8 .357 .932 3.54 
3 12.3 4.2 .238 902 4.06 
4 8.1 1.2 .833 875 .88 
5 6.9 2. 5.000 .863 13 
6 6.7 1 10.000 .64 
7 6.6 ) 0 ) 
Total loss: 10.4 Average: .909 11.87 
(26.0/ha) (29.675/ha) 


Table 2.—con. 


Nes i gy Ed RI en Nn ee a Oe 
Number of trees 


Diameter Year of per acre Predicted 
class infestation Gree aad 1/D, a tree loss 
pene iniekes weer entree On ts A ee et we ee i 

14-inch 0 8.0 0.4 2.500 0.880 0.89 
(35.6-cm) 1 7.6 Wee .833 .867 2.28 
2 6.4 2.2 454 .826 3.07 
3 4.2 2.0 .500 724 1.88 
4 2.2 4 2.500 .606 25 
5 1.8 2 5.000 55 10 
6 1.6 0 0 0 
Total loss: 6.4 Average: 743 8.47 
(16.0/ha) ; (21.175/ha) 
16-inch 0 2 0.3 3.333 0.598 0.624 
(40.6-cm) 1 1.8 a) 3:338 5.45 534 
2 125 7 1.429 .407 B41 
$} 8 3 Sheil) .209 .238 
4 5 2 5.000 .078 105 
5 .S = 10.000 .017 .033 
6 2 0 0 0 
Total loss: 1.9 Average: .309 2:375 
(4.75/ha) (5.925/ha) 
>16-inch 0 2.0 03 3.333 0.582 0.63 
( >40.6-cm) 1 lve 4 2.500 it .67 
2 es} 9 destitit .270 .88 
3 4 si 10.000 .056 05 
4 3 2 5.000 .004 .06 
5 at 0 0 0 0 
6 a 0 0 
Total loss: 1.9 Average: .285 2.29 
(4.75/ha) (5.725/ha) 
Total 0 211.0 3.0 5.2 
1 208.0 5.0 7.7 
2 203.0 14.0 16.2 
3 189.0 27.0 15.6 
4 162.0 4.0 3.0 
5 158.0 2.0 les 
6 156.0 1.0 2 
U 155.0 0 0 
Total loss: 56.0 Average: 153 59.2 
(140.0/ha) (148.0/ha) 
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Figure 1.—Predicted versus observed tree loss by 
year based on q; by 2-inch tree diameter class 
(Situation A, Klein and others 1978). 


Table 3.—Predicted versus observed tree loss by year based on q;, the average probability of tree loss, by tree diameter 
class and stand (Situation A, Klein and others 1978, grouped by greater diameter spread) 


Number of trees 


Diameter Year of per acre Predicted 
class infestation Groen - “Dead 1/D, ot tree loss 
a a eR ee ea ee 
6-12 inch 0 198.3 2:5 0.4000 0.997 25) 
(15.2-30.5-cm) 1 196.8 3:2 3125 .995 3.1 
2 193.6 9.7 .103 995 9.2 
3 183.9 25.0 .040 994 2a 
4 158.9 3.3 .303 .994 2.6 
5 155.6 eS .769 994 1.0 
6 154.3 a) 3330 0 a2) 
7 154.0 0 0 
Total loss: 44.3 Average: .995 40.3 
(110.75/ha) (100.75/ha) 
14-16 inch 0 12.1 1.0 1.000 0.917 2:5 
(35.6-40.6-cm) 1 Ailed 1.9 526 .906 3.9 
2 9.2 3.8 .263 .869 5.4 
3 5.4 2.4 .416 183 2.3 
4 3.0 8 1.250 679 5 
5 22 3 3.333 613 Alb 
6 1.9 0 0 0 0 
Total loss: 10.2 Average: 194 14.8 
(25.5/ha) (37.0/ha) 
>16-inch 0 4.1 0.6 1.667 0.768 16 
('>40.6-cm) 1 3.5 WA 1.429 127 Vat 
2 2.8 1.6 .625 589 2.0 
3 1.2 4 2.50C 363 3 
4 8 4 2.500 177 2 
5 4 4 10.000 .056 03 
6 13 0 0 0 0 
Total loss: 3.8 Average: 447 5.9 
(9.5/ha) (14.75/ha} 
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Figure 2.—Predicted versus ob-served tree loss by 
year based on q;, by grouped tree diameter 
classes and total stand (Situation A, Klein and 
others 1978). 


28 All Trees 


Situation A 


TREE LOSS 


YEAR 


Figure 3.—Predicted versus ob- 
served tree loss by year for total 
stand based on q; (Situation A, 
Klein and others 1978). 


The second data set came from a mountain pine beetle 
infestation in the Gallatin River Drainage, and was used 
only for total tree loss, because the data were not 
originally taken by diameter classes (situation B). Tree 
loss over time did not fall into the usual bell-shaped pat- 
tern, yet the predicted tree loss approximated the actual 
double-peaked curve (table 4; fig. 4) (Burnell 1977). 

Answers to the questions, ‘‘Which of the lodgepole 
pine stands are the most susceptible to mountain pine 
beetle outbreak development?’ and ‘‘How many trees 
will the manager lose if the stand becomes infested?’ de- 
pend upon risk. A definition of risk has two parts: (1) 
probability of an outbreak within a set time period, and 
(2) expected loss in the advent of an outbreak (Safranyik 
1982). Reliable methods are not available to predict 
when an outbreak will develop, but we can predict the 
most susceptible stands and also forecast stand deple- 
tion in terms of stand structure should an epidemic oc- 
cur. To date, six risk-rating systems have been 
developed that are based on climatic and tree/stand 
variables having a major effect on beetle survival and 
distribution. Rate of spread could be considered using 
historical maps (fig. 5) or mathematical models based on 
habitat type. 


Table 4.—Predicted versus observed tree loss by year based on q;, the average probability of tree 
loss by tree diameter class and stand (Situation B, Burnell 1977, grouped by total stand) 


Number of trees 


Year of per acre Predicted 
infestation Green Dead 1/D, at tree loss 
0 370.2 3.6 0.278 0.997 6.6 
1 366.6 8 1.250 .997 1.5 
2 365.8 19.5 0513 .997 34.1 
3 346.3 16.4 .0609 .997 2723 
4 329.9 77.8 .0128 .999 106.5 
5 25244 31.8 .0315 .992 37.1 
6 220.3 10:3 .0971 .987 aaa 
7 210.0 0 0 0 0 
Total loss: 160.2 Average: .995 224.2 
(400.5/ha) (560.5/ha) 
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Figure 4.—Predicted versus ob- 
served tree loss by year based on 
gq, by tree diameter class, grouped 


and total stand (Situation B, 


Burnell 1977). 
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Figure 5.—Chronology of mountain pine beetle infesta- 
tions, Glacier National Park and Blackfeet Indian Reser- 
vation, Mont. 1972-1980. (McGregor and others 1982) 
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Historical maps are useful in areas that have suffered 
repeated severe outbreaks where remaining stands can 
be hazard rated. Through forest inventory and survey 
data, forest cover types can be delineated showing 
mature, overmature, pole, and reproduction-size stands. 
When location and stand composition are known, maps 
can be composed depicting susceptible stands. These 
maps provide a rough hazard rating over large areas, 
which managers can use to initiate strategies to prevent 
future infestations or to salvage logs and reduce fuel 
loads in stands devastated by the mountain pine beetle. 
Usually managers can expect that another epidemic will 
begin within 20 to 40 years, when remaining trees reach 
size classes with phloem thickness conducive to a 
population buildup (Amman 1975). However, this 
depends on characteristics of stands and how soon 
residual trees become susceptible; and it is likely that in- 
festation recurrence will be prolonged in managed 


stands. Historical maps, timber type maps, and timber in- 


ventory surveys can provide the basis for hazard rating 
stands. The ratings can be in very broad, but also ex- 
tremely accurate, categories (McGregor 1982). However, 
significant differences occur within areas rated high 


11 


hazard as to the amount and rate at which mortality 
develops, peaks, and subsides in various stands. 
Relating mortality with habitat type on a stand basis 
has helped refine hazard rating of lodgepole pine stands 
in the Forest Service Northern Region. 

The Insect and Disease Damage Survey Model 
(INDIDS) (Bousfield 1981) and our Rate of Loss Model 
were tested using approximately 1,200 stands with 
mountain pine beetle infestations ranging from | year to 
the end of the epidemic (McGregor and others 1982.) The 
INDIDS Model is used to analyze forest insect and 
disease data collected from variable or fixed plots. It 
uses summaries of detailed mensurational data of in- 
fested and residual green stands—a tree species, size 
class, and damage class for each designated survey type. 
Use of INDIDS Model results in computations of tree 
and volume loss and basal area killed per acre (Dilworth 
and Bell 1968). 

The Rate of Loss Model was integrated with the 
INDIDS Model to estimate mortality trends for infested 
stands or to obtain loss estimates (tree, cubic, and board 
foot volume) by diameter class over infestation time for 
green stands, should they become infested (table-5). 
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Figure 6.—Predicted lodgepole pine trees and volume loss from mountain 


pine beetle by habitat type over time for Madison Ranger District, 
Beaverhead National Forest, and Hebgen Lake Ranger District, Gallatin 


National Forest, Mont. 


18 


Omrammw OA AWOrF MZCrod NO AMM OHDBCO NO x OMramMmMwm COA 4HOr MxrCro< NO AMM OHMWCO NO x* 


onmrsammo O4 4HOr macro< NO AmmMnN OHWECO NO x 


MADISON 
Y=-8 36772. 438K PSME/VASL-VAGL 
Y=-8 .274+8.434"X PSME/SYAL 
Y=-2.285+2 .322"X PSME/SYAL—CARU 
Y=-@_133+23.1SSexX PSME/SYAL—-SYAL 
Y=-@ . 437+8.4S3»X PSME/CARU 
Y=-2 476-2 _S4SeX PSME/CARU-CARU 


MADISON 
Y=—-B.327+8.327"xX PSME/CAGE 
Y=-2.887+8.487"X PICEA/LIBO 
Y=—-B.835+28.S2B"x ABLA/VACA 
Y=-28.785+2.732"xX ABLA/CACA 
Y=—®.668+2.733"xX ABLA/LIBO-LIBO 
Ye-® .482+28.631"X ABLA/LIBO-VASC 


HEBSEN 
Y=—-8.888+28.458"X ABLA/VAST-VASC 
Y=B.S84+2.142"X ABLA/ALSI 
Y=8_184+8.196"X ABLA/CARU 
Y=. O9B2-B.B18eX ABLA/CAGE 
Y=O.823+8.835eX ABLA/CAGE-PSME 


MADISON 
Y=-@.252+08.277"X ABLA/VAGL 
Y=-8.387+8.388"xX ABLA/VASC 
Y=-8.229+0.264"X ABLA/VASC-CARU 
Y=-9.265+8.286xX ABLA/VASC-VASC 
Y=-8.248+9.268"xX ABLA/ALST 


OmMramme® of AHOr umm~na TNO x 


MADISON 
Y--@.185+.209"x ABLA/CARU 
Y=-8.277+8.277"X ABLA/CAGE 
Y=-@.155+0.1558X ABLA/CAGE-PSME 
Y=-0.348+0.354"x ABLA/PIAL-VASC 
Y=-8. 153+. 163~x PICO/CARU 


oOmramm®@ O4 AHOr HMMA TO x 


Figure 6.—(con.) 


Stand data were then subjected to analysis of variance 
and analysis of covariance for completely randomized 
design and graphed to show lodgepole pine mortality by 
habitat type over time (fig. 6). Analysis shows that the 
percentage of lodgepole pine killed and volume loss vary 
by habitat type. 

In some habitat types, tree mortality increased rapidly 
and most susceptible trees and all volume are killed in a 
relatively short time (fig. 6; ABLA/VASC-VASC, 
ABLA/ALSI). In others, mortality may occur over a 
10-year period and never exceed 30 percent of the stand 
(fig. 6, ABLA/CARU, ABLA/LIBO-LIBO). All suscepti- 
ble trees may be killed in other habitat types, but it 
may require 8 to 10 years. These data provide guidance 
as to which stands within those classed as high hazard 
should receive priority management. For example, 
management may be postponed until the next decade if 
stand mortality does not exceed 20 to 30 percent over a 
10-year period. Meanwhile, stands can be rated and 
management implemented in the stands containing 
habitat types where considerable tree mortality or 
volume loss is predicted to occur over a short time. By 
putting the higher risk stands under management, loss 
would probably be prevented in some high-, many 
moderate-, and many low-risk stands. 
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INTEGRATION WITH FORPLAN 


The Forest Service currently uses FORPLAN, a linear 
programing model (Johnson and others 1980), for land 
management planning which is the basis for land use 
allocations and scheduling of management activities. The 
management activities and associated products, costs, 
and environmental effects used in FORPLAN are 
reflected in prescriptions for stands within analysis 
areas. In the Forest Service Northern Region, analysis 
areas are lands that meet certain common classification 
criteria; these lands are not usually contiguous. 
Classification criteria include habitat type, timber size 
class, slope class, and other characteristics. Prescriptions 
describe specific management practices used to manage 
specific stands. 

One approach to modeling tree mortality caused by 
the mountain pine beetle using FORPLAN has been to 
predict susceptible areas within analysis areas and prob- 
able mortality over two decades. Although it might be 
possible to predict rate of loss caused by the beetles 
throughout a forest, this information would be of little 
value for adjusting yield tables if the locations of high-, 
moderate-, and low-risk stands are not identified within 
analysis areas. The FORPLAN model would spread bark 


beetle effects over the next two decades for all stands 
within analysis areas, which would not allow scheduling 
earlier harvest of stands with a high probability of in- 
festation and mortality. 

Another approach is recommended when the location 
of stands within analysis areas is identified through 
timber surveys or stand examinations. Beetle attack 
may then be simulated by a “‘prescription”’ that shows 
the effects of an epidemic in the absence of timber 
management. If other management practices were not 
implemented, it would be necessary to constrain the bee- 
tle “‘prescription’”’ by assignment to a certain acreage. 
Thus there would be two prescriptions—one for some 
stands in parts of the analysis area with infestation, and 
one for other parts with no effects of infestation. 

Stands in the Helena National Forest were analyzed in 
a FORPLAN run by grouping habitat types so mortality 
factors could be directly applied to yield tables. A pro- 


cedure was adopted and used to adjust yield tables based 


on the coefficients developed for the Helena National 
Forest plan (Brohman and others 1982). Coefficients 
were based on the assumption that a 50 percent loss of 
lodgepole pine would occur over a 5-year period. The 
estimated loss as a percentage of volume by age classes 
was determined as shown: ; 

Miao tl cL) 

Y,"= Y,(1 — % L) 

Nis Non (e le) ey ceed 


where 


= proportion of volume lost to beetles (50 percent = 0.50), 


Y, = tabular volume for decade j of the plan, and 
Yj;' = adjusted volume expected to exist in decade j. 


Such coefficients must be derived for each habitat 
type or habitat type group to be applicable to the model. 
Decade 1, 2, or 3 of the Forest Plan may correspond to 
different decades in the yield table for different stands 
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or habitat type groups within analysis areas. For exam- 
ple, if groups of stands are 105 years old, then Y, is the 
tabular yield shown at 110 years (25 percent loss by 
year 5). If the current age is 165 years, then Y, is the 
tabular yield shown for 170 years (25 percent loss by 
year 5, and 75 percent loss by year 10 at 175 years). The 
graphs in figure 7 were developed using this approach 
and the INDIDS/Rate of Loss Model for the Helena Na- 
tional Forest in the absence of beetle attack. The factor 
or proportionality is (1—L), the proportion of stand 
volume not killed. 

That the predicted results graphed in figure 7 will ac- 
tually happen is questionable. Beetle-induced mortality 
will reduce competition for trees that are not attacked. 
However, trees not susceptible to bark beetle attack are 
usually smaller and less vigorous. These trees will pro- 
bably respond to a decrease in competition. But amount 
of response will depend on tree age and various site fac- 
tors. We do not know at what rate the remaining live 
stand will grow compared to what it would have done 
without attack. 

The final step in the FORPLAN run for the Helena 
National Forest plan was to adjust existing yield tables 
by the appropriate coefficient for each habitat type 
group. Regenerated stand tables were not adjusted, 
because management should prevent mountain pine bee- 
tle outbreaks over a rotation. The assumption that the 
beetle will infest susceptible stands throughout the 
National Forest in the next 20 years may not be totally 
correct, but it seems probable based on available infor- 
mation. By including coefficients in the yield tables, the . 
FORPLAN model should show which highly susceptible 
lodgepole stands need immediate harvesting, and which 
stands should be harvested before becoming highly 
susceptible. By using assessments from FORPLAN and 
harvesting in high-hazard, susceptible stands before an 
epidemic develops, land managers should be able to 
minimize tree mortality. caused by the beetle. 


DRY FIR FOREST 1 r DRY FIR FOREST 
Cler/ T My or Yo-8.4120+2,49884xX DF/CARU 
_————— Ys-8.2312+8,2680"xX DF/CARU She Ye-@. 1438+2.28864X DF/PHMA 
oa eS eS Y=-8.9452+8.1228"X DF/PHMA OU aoa res ate ane Yo-@8.2481+8.4495"xX AF/CARU 
eee er eee Ye-@. 1S68+8.2348"X AF/CARU ee oe 
8.6 c 
% If u 
B 8.8L 
I 
2 c 
7 oS bt (te tcine ee 
R E 
E E 
3 T 
s 8. 
a4 ‘ 
L 
ry re 
Ss Ul, «ay 
t M 
i= 
T 
L. 42! 
2 oO 
Ss 
E T 
E a. 
Th T 
& ie) 
5 B a. 
i= 
— 
T 
ee Ss 
E 
Ss 
2 4 =e 4 1 = 4 aa 4 4 
@ 1 cA 3 4 Ss 6 Pe 8 9 18 
YEAR 
a 
(haf hs COOL SLOPES 1 r COOL SLOPE 
= % ——_- Yo-8.3737+8.5185"xX AF/LIBO 
———— Y=-8.1883+8,2SS2@xX AF/LIBO ee ee Ye-9.2476+8.4716"X AF/VAGL 
— i Yo-@8.1379+8.2218xX AF/VAGL fe} Yo-@8.5363+2.6391X AF/VASC 
ty. ae Y=-@,2655+8,2978"X AF/VASC Fo agL Yo-94878+0_6169NX AF/MEFE 
@.6/L ¥=-8.3656+0.4181~x AF/MEFE c 
x” U 
B 8.6L 
a I 
F c 
sy ghsasey Fi <8; 
R E 
E E 
(= a9 
s 8. 
8. Vv 
L oO 
oO L 
Ss US8. 
T M 
E 
etic 
fe kL. 8. 
Oo 
B s 
E ag 
E 8. 
7 38. t 
L oO 
E 
s B @. 
E 
e. E 
T 
iE. 4a. 
E 
s 
@ a = i i 4 4 4 4 4 4 —t _! 
@ 1 2 3 4 Ss 6 7 8 3 18 
YEAR YEAR 
8.7 r FOREST TYPE % 1 y. ‘4 ene ORY eee RAL 
a 4 =-B8.26414+8.4313mx R FORES 
es Y=-8.1528+8.2228"xX DRY FIR FOREST Y=-8_3793+8.54@4"X COOL SLOPES 
= She Y=-@.2838+8.26S6"xX COOL SLOPES le) Ys-@8.2482+8.4199"xX MESIS S-S FOREST 
ae Yo-8.1494+8.2662"x MESIS S-S FOREST Fo a. 
8.6 L c 
4 U 
B a.eL 
oO xi 
F c 
Tames F 8.7 
R E 
E E 
E T 
s 8.6 
8.4 v 
{i oO 
oO (ls 
s U 3@.sS 
A M 
E 
Gees 
fe) Ll @.4 
o 
B s 
E T 
E 8.3 
Tae T 
L o 
(= 
s B 8.2 
E 
a. E 
T 
a) | 
cael 
s 
a L — ee 4 4 4 —t 1 — 4 =) a 4 1 ie =i rt 1 4 4 1 = 
8 1 2 3 4 5 6 oh 8 3 18 a8 1 2 3 4 s 6 7 8 3 12 
YEAR YEAR 


Figure 7.—Projected tree and volume loss from mountain pine beetle for 
lodgepole pine habitat type groups within dry fir, cool slopes, and mesic 
sites on the Helena National Forest. 
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Forest managers have less than a 50 percent chance of growing lodgepole 
pine to 16-inch diameters in unmanaged stands because of recurrent depreda- 
tion from the mountain pine beetle. Hazard rating methods provide techniques 
for managers to identify susceptible stands. The Rate of Loss Model refines ex- 
isting risk rating systems, and provides a method for predicting tree and voiume 
loss by habitat types. This model is provided to assist land managers in projec- 
ting tree mortality over time, and as a link with the FORPLAN Model for use in 


forest planning. 
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